Dissociation of Vortex Stacks into Fractional-Flux Vortices 
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We discuss the zero field superconducting phase transition in a finite system of magnetically 
coupled superconducting layers. Transverse screening is modified by the presence of other layers 
resulting in topological excitations with fractional flux. Vortex stacks trapping a full flux and present 
at any finite temperature undergo a dissociation transition which corresponds to the depairing 
of fractional-flux vortices in individual layers. We propose an experiment with a bi-layer system 
allowing to identify the dissociation of bound vortex molecules. 



The zero field superconducting to normal transition in 
thin films and layered superconductors is triggered by the 
proliferation of topological defects; the unbindingof Pearl 
vortices Q in thin films and of pancake vortices [3 in lay- 
ered superconductors generates a Berezinskii-Kosterlitz- 
Thouless transition (BKT) 0] which has been studied in 
detail 0,0] ■ New features emerge when going to a system 
with a finite number N of magnetically coupled layers. 
Besides Pearl type vortex stacks penetrating through the 
full array of layers, cf. Fig. 1(a), fractional- flux vortices 
appear which reside within the individual layers 0, , 
the analogue of the pancake vortices in a layered material; 
the reduced trapped flux associated with these vortices is 
due to the presence of other layers modifying transverse 
screening in the multi-layer system. The Pearl vortex can 
be viewed as a linear arrangement (stack) of fractional- 
flux vortices; the intra-layer unbinding transition of these 
fractional-flux vortices then corresponds to the dissocia- 
tion of full-flux vortex stacks present at any temperature, 
cf. Fig. 1(b). In this letter, we discuss the prospects of 
observing this dissociation transition in an experiment; 
in particular, we study a bi-layer system in a counter- 
flow geometry which allows to observe the dissociation of 
vortex molecules into half- flux vortices, cf. Fig. 1(c). 

The basic prerequisite for the appearance of a BKT 
transition is the logarithmic interaction between defects, 
V{R) = 2e2 In R, where we attribute an effective 'charge' 
e to the defects. In the absence of screening, e.g., in 
a superfluid ^He film, the logarithmic interaction be- 
tween vortices extends to infinity and the system un- 
dergoes a finite temperature BKT transition In a 
superconducting film, transverse screening restricts the 
log-interaction to the screening length A; the Pearl vor- 
tices assume a finite self-energy V{A)/2 and hence can 
be thermally excited at any finite temperature — the 
superconducting to normal transition then is shifted to 
T = 0, although a sharp crossover survives at a finite 
temperature T^^,.^ « e^/2 0. Below, we concentrate on 
a system with N magnetically coupled superconducting 
layers, i.e., vortices interact through the transverse mag- 
netic potential A, while the Josephson coupling due to 
Cooper pair tunneling between the layers is assumed to 
be negligible, as it is the case in a material with insu- 
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FIG. 1: (a) Fractional flux vortex and vortex stack in a A'^ = 5- 
layer system, (b) The excursion of one ff vortex from the vor- 
tex stack is equivalent to the combination of a complete stack 
and a vort ex-ant i- vortex pair; pair unbinding in a layer is 
equivalent to stack dissociation, (c) Geometry for the bi-layer 
setup shortcircuiting the effect of vortex stacks, (d) Sketch of 
I-V curves at various temperatures (see arrow); the algebraic 
dependence 1/ cx with a > 3 at low temperatures turns 
ohmic [a = 1) at high temperatures. The regime between 
Tbkt S'lid Tbkt contains the interesting features associated 
with half-fiux vortices. The current scale loft ~ -^o(C/-^cff) 
(vertical dashed line) separates the physics of unscreened vor- 
tices from that of half-fiux vortices. At T^k^ the I-V curve 
exhibits the characteristic exponents 3 and 5 at small and 
large currents tracing the crossover from unscreened to half- 
fiux vortices. Above Tbkt additional ohmic regime due to 
free half-fiux vortices appears at low currents; the dotted line 
marks the crossover current 7hf- At temperatures T > Tbkt 
the ohmic regime takes over and leaves only a small non-linear 
region at high currents probing unscreened vortices. 



lating layers separating the superconducting ones. The 
presence of additional layers leads to drastic modifica- 
tions in the potential Viji) between individual vortices in 
the same layer: i) the log-interaction extends to infinity, 
ii) the magnetic flux <I>t trapped by individual vortices is 
reduced to a fraction 1/A^ of the flux unit $o = hc/2e, 
$t(Af) = ^q/N; the effective coupling in the inter- 
action potential between defects is reduced correspond- 
ingly, ^ e^[l — 1/-/V]. Hence, the additional layers act 
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on the vortices in the same way as a dielectric matrix acts 
on charged particles. The unbinding of the fractional-flux 
(ff) vortices triggers a finite temperature BKT transition 
at Tbkt ~ e^[l— 1/-/VJ/2. At the same time, vortex stacks 
are present at any temperature; the unbinding of (ff ) vor- 
tices at Tbkt then describes the dissociation of the vortex 
stacks rather than the superconductor-normal transition 
which, strictly speaking, appears already at T = 0. Fi- 
nally, in a bulk layered superconductor, the magnetic 
field escapes in the transverse direction and the flux 
trapped by pancake vortices saturates at $t = rf'&o/SA, 
where d and A denote the layer separation and the bulk 
planar penetration depth. The pancake-vortex unbind- 
ing at Tbkt ~ e^/2 describes a generic phase transition 
as no competing stacks show up in the bulk. 

Below, we address the prospect of observing fractional- 
flux vortices in an experiment. The presence of vor- 
tices and their unpairing reveals itself in the current- 
voltage characteristic. Fractional-flux vortices appear 
most prominently in a bi-layer system, see Fig. 1 (c) . We 
propose a counterflow experiment where the contribution 
of stacks (vortex molecules) is eliminated. The current- 
induced unpairing of vortices then produces an algebraic 
characteristic V oc I"; the change in the exponent from 
a = 3 to 5 with increasing current, cf. Fig. 1(d), is a man- 
ifestation of the stack dissociation transition at T^kt- In 
the following, we briefly derive the structure of topolog- 
ical excitations in layered systems, discuss their thermo- 
dynamics, and analyze the features in the I-V character- 
istic related to fractional-flux vortices. 

Consider a superconductor of thickness d and with a 
London penetration depth A. Central to our discussion 
is the interaction potential V{R) between vortices: The 
current associated with a vortex is driven by the 27r-phase 
twist Vip = -n,xR/R'^,j{R) = -(c/47rA2)[V<y9$o/27r + 
A(i?)]; transverse screening through the vector potential 
A reduces the action of Vip until complete compensation 
is reached once a full flux quantum <i>o is trapped. A 
second vortex placed a distance R away is subject to the 
Lorentz force F{R) = -j^{R)d^o/c = -{2eod/R)[l - 
$(i?)/<I>o], where $(i?) = 2-kRA^(R) is the flux accu- 
mulated within the distance R and Eq = (4'o/47rA)^ is 
the line energy; integrating this force provides us with 
the desired potential between the vortices. The incom- 
plete asymptotic screening with a reduced trapped flux 
4>t = $(i? oo) < $0 then gives rise to a logarithmic 
interaction ViK) ^ 2eQd\V — ^t/'&o] ln(i?) and hence a 
BKT phase transition (note that ^ Eqc')- 

In order to flnd the flux <I>t we have to solve the 
Maxwell equations for the potential A. We consider a 
stack of N superconducting layers of thickness d^ and 
separated by a distance d. The penetration depth As 
of the layer material deflnes the bulk planar penetration 
depth A^ = Agd/dg. We place the vortex at the origin of 
the layer positioned at z = and describe the protecting 
layers of thickness c?^ = (J>S — n)d and d^ = {n — l)d 



above and below the fllm in a continuum approximation. 
The system 

V2A-^A=^(A+gv^)j(z), -d<<z<d> 

V^A = 0, z < -d" and z > d> (1) 

then assumes the solution 

^od rdKMKR),,^ . 
A,iR,z)^^J^ ^^^fiK,z), (2) 

with Ji the Bessel function and the function f{K, z) — 
[1 - a<i'(iir)]e"-^+l"^l + ad'{K)e^+^''^ describing the z- 
dependence within the superconductor. Here, ad' (K) = 
( K+ - K)/ [{K+ + K)e^^+'^' + {K+ - K)], with K+ = 
y/K'^ + A^2 and d' = d< {d' = d>) in the region < 
z < (0 < z < d^). The denominator C{K) assumes the 
formC(i^) = [l-2ad<{K)]K+ + [l-2ad>{K)]K++d/X^. 
Outside the superconductor, the fleld is obtained by re- 
placing z > (z < 0) by {~'d,'^) in / and an additional 
correction factor exp[K{d^ — z)] {exp[K{d^ + z)]). The 
Pearl- and pancake vortices are recovered in the limits 
d<, d> =0 and d<, d> = oo. 

The magnetic flux $t trapped by a vortex is extracted 
from the vector potential at z = 0; for a thin iV-layer 
system with Nd <C A we find the asymptotic form 

A^R > Acff, z = 0) - ($t/27ri?)(l - X,s/R) (3) 

with parameters Acff ~ 2\^/dN and <i>t(A^) ~ ^o/N 
independent of the layer position n. For large N the 
trapped flux saturates at the value $t = d^Q/2\ and 
Acff = 8], while $t = and AoS = A = 2X1 /ds 1°^' 
thin film. The decrease $t = 'I'o hi trapped flux with 
increasing number of layers is easily understood: with 
Nd <C A, the same magnetic field (and hence the same 
fiux) is penetrating the N layers. On the other hand, a 
vortex stack (i.e., N (ff) vortices) carries a full flux $0. 
The flux associated with one individual (ff) vortex then 
is a fraction 1/iV of the value trapped by the vortex stack 
and thus (f>t = ^o/N. 

The incomplete screening of the vortex singularity 
produces a log-interaction between vortices, V{R) — 
2eod\n[R/£) at small distances R <C Acff and V{R) — 
2eod[ln(Acff/0 + '^t/'^o] ln(i?/Acff)] involving a large 
self-energy but a reduced prefactor at large distances 
R 3> Acff. This logarithmic interaction competes with 
the entropy of vortex-pair excitations and t rigg ers a 
Berezinskii-Kosterlitz-Thouless transition at 0, E3 

In a thin fllm, <I>t = $0 and the flnite range of the inter- 
action between Pearl vortices pushes the real transition 
to T = 0, in agreement with Q; the presence of one ad- 
ditional 'protecting' layer changes the situation: trans- 
verse screening reduces the trapped flux to half its value. 
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$t = ^"0/2, thus extending the range of the logarithmic 
interaction to infinity and pushing the transition tem- 
perature to a finite value T^kt = e-Qd/A. Adding more 
(A'' — 1) layers, the trapped flux $t = $0/-^ decreases 
further until assuming the asymptotic value $t =d$o/2A 
in a bulk superconductor where T^kt ~ £0*^/2 is largest. 

The appearance of a finite temperature phase transi- 
tion due to the protecting action of additional supercon- 
ducting layers has its counterpart in multigap supercon- 
ductors 101, cf. also Ref. ^3 In both cases, the creation 
of a topological defect in one superconducting component 
or layer induces screening currents in the other compo- 
nents/layers via coupling to the common gauge field A; 
the resulting incomplete screening extends the interac- 
tion between defects to infinity, although with a reduced 
'charge'. A finite Josephson coupling between the layers 
or between the components of a multigap superconduc- 
tor spoils this phenomenon through the appearance of 
a linear confining potential. While this coupling can be 
(made) arbitrarily small in a layered system, the internal 
Josephson effect in a multi-component superconductor 
cannot be tuned and is not necessarily small Q|. 

The setup where fractional-flux vortices make their 
most prominent appearance is the bi-layer system. Its 
thermodynamic properties are obtained from an exten- 
sion of the renormalization group analysis in a film [fl[llO|| 
and involves the flow of the superfluid density K{R) at 
scale R with K{£) = Kq — eod/irT and the vortex fugac- 
ity y{R), with {y/R)'^ the density of vortex-anti- vortex 
pairs of size R. Here, the renormalization involves a 
two-stage process: i) unscreened vortex pairs are inte- 
grated out on scales R < AcS and provide renormal- 
ized values {K,y) aX R — Acff; ii) the flow is restarted 
with a reduced coupling eod/2 and half-flux (hf) vor- 
tex pairs are integrated out on scales R > Xctf- We 
obtain the following results (see for details): At 
temperatures T < T^kt = £od/4: the initial fugacity 
2/0 = exp(— E'c/T) due to the core energy Ec ~ £0^ 
flows to zero and the superfluid density K Ki K > A/tt 
remains finite. Above T^kt (hf) vortices start unbind- 
ing: a narrow critical regime (with a correlation length 
6if ~ Aeff exp(^/2V6A<), At = (T - Ti'i^)/Ti'i^ and b 
a dimensionless parameter) is followed by a mean-field 
behavior where the fugacity diverges to infinity and the 
renormalized superfluid density vanishes beyond the cor- 
relation length ^hf, defining a density of free (hf) vortices 

with a of order unity. On approaching T^^^ = £od/2 the 
correlation length becomes comparable to Acff; beyond 
'^BKT unscreened vortices start unbinding at small scales 
below Acff. Note that vortex stacks do preempt the super- 
conducting transition of the bi-layer system but preserve 
superconductivity within the individual layers; indeed, 



the force of a (hf) vortex acting on a vortex stack van- 
ishes rapidly beyond the effective screening length Aofr 
due to the complete screening of vortex stacks. 

The presence of (hf) vortices can be traced in an ex- 
periment measuring the I-V characteristic (we denote 
the sheet-current density by / = jd). In the counter- 
flow geometry of Fig. 1(c), the applied dc current acts 
oppositely on the two (hf) vortices constituting a stack 
and the linear response due to drag motion is quenched. 
The current-induced dissociation of pairs and stacks of 
(hf) vortices produces a non-linear I-V characteristic; the 
change in slope from 3 at low currents to 5 at high cur- 
rents signals the thermodynamic dissociation of stacks 
at Tbkt- The current-induced unbinding of (hf) vor- 
tex pairs involves a thermal activation over the bar- 
rier U{I) = uia.xR[V{R) - I'PoR/c] « 2£orfln(i?//0 at 
small distances Ri <^ Acff, while U ~ 2eodln(Acff/C) + 
eod\n{Ri / Xos) for Rj » Acff. Here, Ri ~ S,Iq/I de- 
notes the unbinding scale and Iq = 2eodc/^o£, is close 
to the depairing current. Applied currents smaller than 
/off = /oC/AefT probe lengths larger than Aefi and the 
effects of half-flux vortices become accessible. 

The equilibrium density of free vortices derives 
from the steady state solution of the rate equation 
dtUy = r — ^^n^/Trcc, with F oc exp(— [//T) the produc- 
tion rate of free vortices and ^'^ /t^cc the recombination 
parameter. Vortex drag then produces a finite non-linear 
resistivity, p « C^Pn^^v with the normal state resis- 
tivity, and a corresponding algebraic I-V characteristic 
V/Vo - (///o)"('^' with 

a(r,/) = l + 7rX(T)[l-$t(i?/)/$o]. (6) 

The exponent a depends explicitly on the flux associated 
with the vortices: at short scales, unscreened vortices are 
probed and a = 1 + nK. On the other hand, large dis- 
tances probe half-flux vortices and the exponent is re- 
duced to a = 1 -I- ttK /2. The crossover between these 
two regimes appears at the current /off ^ -^o- This re- 
duction in a at /cff is the most prominent feature in the 
I-V characteristic signalling the presence of half- flux vor- 
tices in the system; at T^kt) the change in slope is from 
5 at large currents to 3 at low currents, cf. Fig. 1(d). 
The associated voltage signals are weak as the density of 
free vortices is already small, ^ l/A^fj at T^^^ and 
nv ^ ^^/Agfj at the true transition point Tbkt- 

The temperature T^kt defines a resistive transition 
due to the proliferation of free half-flux vortices. Above 
Tbkt the ohmic resistance appearing at low currents 
I < Im — Io£,/£,hi (probing distances larger than the 
correlation length ^hf) is determined by the density of 
free (hf) vortices, p « PnC^'^hf ■ In the mean-field regime 
above Tbkt we can make use of (0 and find the loca- 
tion of the crossover at /hf ~ Io{Xcif / i){i^/yo/ XcsY^'^^*' 
and the temperature dependent resistivity pm{T) « 
/On(Aoff/C)^(?\/yo/'^cff)°^^* due to free half-flux vortices. 
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The measurement of /hf or phf in this regime provides di- 
rect access to the correlation length ^hf and its mean-field 
like temperature dependence. 

In order to analyze the BKT transition at T^kt var- 
ious experimental constraints have to be accounted for: 
i) The system must be larger than the screening length 
Aoff beyond which (hf) vortices appear. The relation 
(El) impHes that XesiTi^^T) ~ 0.5 cm/(Ti'i^T in K), hence 
Aeff « 1 mm in a typical low Tc material, ii) Both 
the interlayer distance d and the layer thickness dg have 
to be small compared to the bulk penetration depth A. 
in) The Josephson coupling has to be small enough to 
push the confinement length — \/ {jo/ 3j)^d beyond 
Aeff; separating the superconducting layers with an in- 
sulator the Josephson current jj can be made ar- 
bitrarily small, iv) The mean-field temperature depen- 
dence oc (1 — T^/r^^)~^/^ of the parameters A and ^ tends 
to push the temperatures Tbkt and T^j^^ towards the 
mean-field critical temperature {T^kt ^ Tbkt)/?c ~ 
(Te - T^^^)/T, « 4«<^°' with Gz(^°' = = 
0)d ^ 1 the two-dimensional Ginzburg number [l^- A 
large Ginzburg number helps in distinguishing between 
the temperatures where (hf ) and unscreened vortices un- 
bind, v) The features in the I-V characteristic identify- 
ing the presence of (hf) vortices involve vortex densities 
which are suppressed by the small parameter ^/AeS . Cor- 
recting parameters for the intrinsic dirtiness of thin films 
(see |l9|, we assume a mean free path I limited by the 
layer thickness ds) we obtain the estimates 



-^BKT 



1.5 • 10-4 T,i/2d3/2/Aco, (7) 
« 3.2- 10-9 TeA^Vco/di (8) 



where all lengths are measured in A and temperatures 
in Kelvin. The results Q and ^ tell us that given the 
(clean-) material parameters Aco and ^co it is not possible 
to maximize both Gi'^°* and ^/Aeff simultaneously by 
varying the thickness ds- A reasonable compromise can 
be achieved if we choose a material with Aco ~ Cco ^ 
1000 A, Tc ^ 10 K, and a thickness ds ^ 500 A; this 
yields ^/Aes ~ IQ-^-^ and Gi<'°' ~ IQ-'^. The small 
value of ^/ Aeff implies a small vortex density and requires 
a high voltage resolution, while the smallness of Gi'-™^ 
requires a temperature resolution in the mK range. The 
characteristic halving in a — 1 signalling the presence of 
(hf) vortices below Jeff involves voltages with log(y/Vb) 
between -7.5 and -12.5. With Vq PnjoL ~ 10 mV 
(we assume pn ^ 100 pQcm, jo(?BKT) ^ 10^ A/cm^, 
and L ~ 1 cm) we find that an experimental voltage 
resolution in the sub-pico-Volt regime [23| allows to trace 
this halving in a — 1 over a substantial temperature range 
below Tbkti although the observation of this crossover at 
Tbkt itself pushes the limits of present days experimental 
capabilities. Alternatively, one can trace the presence of 
(hf) vortices by measuring the characteristic mean-field 
type resistivity Phf(T) below T^-^^^ and through direct 



observation with a scanning SQUID microscope [2l|. 

An interesting analogy appears when comparing the 
present system with the bi-layer quantum Hall setup at 
total filling 1^=1. The latter is expected to undergo a 
BKT transition into an interlayer phase coherent state, 
even in the absence of any tunneling between the layers 
[2^ . The (hf ) vortices discussed above (existing in one of 
two layers and with ± vorticity) correspond to topologi- 
cal excitations (merons) with charge ±e/2 and ± vortic- 
ity ,23j ; bound neutral meron pairs have their analogue in 
intralayer vortex-anti-vortex pairs, while bound charged 
merons correspond to vortex stacks. The unbinding of 
meron-pairs in the BKT transition destroys the interlayer 
phase coherence and can be traced in the same type of 
counterflow experiment |24l| as discussed above. 
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